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The model of particle production by strong external sources 
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Abstract 

Using some knowledge of multiplicity disributions for high energy reactions, it is possi- 
ble to propose a simple analytical model of particle production by strong external sources. 
The model describes qualitatively most peculiar properties of the distributions. The gen- 
erating function of the distribution varies so drastically as it can happen at phase transi- 
tions. 
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1 Introduction 



Multiparticle production is the main outcome of high energy collisions. Quantum chromody- 
namics (QCD) provides the general framework for its description. It is especially successful in 



predicting the properties of hard processes. The perturbative approach with some additional 
phenomenological assumptions is effective there (for reviews see [1, 2, 3]). The non-perturbative 
methods are waited for because the bulk of produced partons have a small fraction x of the 
longitudinal momentum of the incoming partners. It was proposed [4, 5, 6] to describe them 
as a classical color field rather than as particles. The large x partons are then considered as 
color sources for the classical field. Many QCD studies of this idea have been done (for recent 
^ \ review see [7]). 

To avoid technical complications of gauge theories, it is proposed [8, 9] to consider the 
simpler theory of a real scalar field <p with strong time-dependent external sources j. The 
^ . Lagrangian density of the theory is 



C = \d^4> - \m 2 <f - |0 3 + j<f>. (1) 



It represents the toy model of the Color Glass Condensate formalism. One hopes that its 
qualitative features could also be valid for particle production in QCD. Moreover, one can use 
the special behavior of multiplicity distributions, gained from QCD studies and confirmed by 
experiment, to develop the more detailed model in the framework of such a theory. This is the 
main goal of the present paper. 



2 General relations 

The multiplicity distribution is the most general feature of multiparticle production processes. 
Any other inclusive (e.g., rapidity) distribution is obtained by averaging over different multi- 
plicities. To be successful, any phenomenological model of these processes should, first of all, 
fit the experimental values for probabilities of n-particle production P n or, equivalently, the 
moments of this distribution. The latter ones are easily computed with the generating function 

oo 

G(u) = J2 Pnu n . (2) 

n=0 
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The unnormalized factorial moments are defined as 



Tq = = £ n{n - i) - {n - q + i)p -' (3) 



n=0 

and the unnormalized cumulant moments are 



d q In G(u) 

}C q = (4) 
They are connected by the iterative relations 

9-1 (<?-!)! 
^ m\{q-m-l)\ 

The mean of the distribution of multiplicities is 

(n) = 5>P n = .Fi = /Ci. (6) 

Factorial moments are always positive. Let us remind that for the Poisson distribution the 
normalized factorial moments F q = J r q /{n) q are identically equal to 1, and the normalized 
cumulant moments K q = IC q /{n) q are equal to except of K 1 — 1. For NBD, both factorial 
and cumulant moments are positive. 

It has been shown [8, 9] that in the theory defined by the Lagrangian density (1) the 
probability to produce n particles to all orders in the coupling g reads 

ft = •-<*'£ = £ P) 

p=l ai+...+a p =n " 

Here a = Y^=iKj P i s the number of disconnected Feynman subdiagrams of the theory pro- 
ducing the n particles, and b r /g 2 denotes the contribution of all r-particle cuts through the 
connected vacuum- vacuum diagrams. 
The generating function reads 

1 oo 

G(u)=exp[-5 5>(« r -l)l- (8) 

9 r =l 

The average multiplicity is 

J oo 

( n > = ^E rfe r- ( 9 ) 

y r =l 

The cumulant moments are 

oo 

^ = - E r ( r - " 9 + !)6r- ( 10 ) 

5 1 r=l 

Cumulant moments of P n play the role of factorial moments of b r /g (compare (10) and (3)). 
The first non-zero terms of the sum begin with b q . 

The "cut" procedure analogous to the Cutkosky rules [10] has been defined [8] for the 
diagrams of the theory with Lagrangian (1). The probability of having p cut subdiagrams is 
equal to 

= (») 

The average number of cut subgraphs is determined as 

(n cut ) = — r (12) 
9 2 
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3 The model 



In principle, the values of b r can be computed [8, 9] with the help of Schwinger-Keldysh for- 
malism [11, 12] or, at leading order, in terms of a pair of solutions of the classical equations of 
motion. At present, this has not been done yet explicitly even though all general expressions 
are given in [8, 9]. Surely, some simplifications will be necessary to solve this formidable field 
theory problem. 

In this situation, one can advocate the phenomenological approach and try to guess such 
behavior of b r which would mimic the shape of the multiplicity distribution favored by QCD 
and experimental data. If successful, this attempt can indicate general features of b r and help 
in approximate computing them from the set of Feynman diagrams. 

The most astonishing regularity of behavior of multiplicity distributions predicted by QCD 
[13, 14, 15] and confirmed by experiment (see review [3]) is the oscillations of the cumulant 
moments as functions of their rank q. According to (10), this peculiar shape can be only 
obtained if b r change their sign with r. To get negative values of b r , the phases of cut diagrams 
in Schwinger-Keldysh approach should play important role. At least for some r, b r can not be 
interpreted as an "intrinsic" r-particle production probability as done in [8]. Correspondingly, 
the relation to the reggeon calculus and AGK cancellations becomes not obvious. At the same 
time, this oscillation should not lead to unphysical negative values of (n) in (9). 

The two-parameter (x, 0) model with oscillating b r , we propose to use, is 



b r = g 2 —rCos(f){r - 1), (13) 



r 

where 



x = -| = P ie a/9 > 0. (14) 



b 
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The generating function and its logarithm read 

G(u) = exp[e micos *cos(xwsin0-0) - e xcos< ^ cos(:r sin - 0)], (15) 

In G{u) = e racos0 cos(a;wsm0-0) - e xcos< ^ cos(:r sin - 0). (16) 

Note the extremely strong and rapidly varying dependence of G(u) on u. Its exponent is an 
exponential function itself multiplied by the oscillating factor! 
The cumulant moments are given by 

K, q = xV^ os ^cos(xsin0+ (g- 1)0). (17) 

They increase exponentially with q and oscillate equidistantly in the g-plane 2 . Their zeros q$ 
are positioned at 

„(Q)_1 | 7r(2m-l) zsin0 

Qm ~ 1 + 20 (18) 

(m is an integer number). The distance between the neighboring zeros is defined only by the 
parameter 

Ag(°) = \. (19) 



2 For NBD, they increase even faster (as q\) and do not oscillate being always positive. 
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The average multiplicity is 

(n) = /Ci = xe XCOB *cos(a;sin0). (20) 

The necessary requirement of positive (n) asks for x sin to be positioned in the hemispheres 
where cos(x sin 0) is positive. Moreover, from the requirement > 1 one obtains for > 

< xsin0 < vr/2. (21) 

The probability of having p cut sub diagrams (11) is determined by 



a 
9 



_ e XCM *cos(a;sin0-0)-cos0. (22) 



Since both (n) and q$ depend on energy, the parameters x and are also functions of energy. 
However, the product x sin must be constant because any its dependence would induce quite 
strange "quasi-oscillating" component in energy behavior of the mean multiplicity. Both (n) 
and qffl should increase with energy tending asymptotically to infinity. For qffl, this property 
follows from QCD asymptotics of K q where the leading order results survive and all cumulant 
moments become positive. Then asymptotically tends to and x tends to infinity so that 
the product X(f) stays (or tends to) constant. For (n) and qf^ one obtains 

(n) = xe x cos(x(p), (23) 

q ? ] = 1 + ^P-. (24) 

If x(f) — > 0, then the asymptotics of (n) and q^ are defined by the parameters x and corre- 
spondingly and, therefore, completely independent. 

The Q CD-like results for the energy dependence of (n) can be obtained with x oc n s 
at s — > oo (up to NLO corrections). For x oc Ins, one gets the power increase of the mean 
multiplicity with energy typical for hydrodynamics [16] and fixed coupling QCD [17, 18]. 

Insering ansatz (13) in Eq. (7) one gets the general form of the multiplicity distribution in 
terms of x and 

Pn = -Jn{<t>), (25) 

where 

n-l 

/n(0) = E c ^ cos ^- (26) 

fc=0 

Unfortunately, no analytical expression for the coefficients c n ^ is obtained. For lowest multiplic- 
ities, they are c 2 , = 1, c 2 ,i = 1; c 3i0 = 0, c 3j i = 3, c 3i2 = 2; c 4 , = 2, c 4 ,i = 1, c 4 , 2 = 11, c 4 , 3 = 4. 

The model is completely fixed if its parameters x and are replaced by the average mul- 
tiplicity and the position of the first zero of cumulant moments. By the order of magnitude, 
these values are similar in different reactions at high energies [19, 20]. To be more definitive, let 
us use our experience in e + e~-collisions at LEP energies [3]. There, the first zero of cumulant 
moments q[ ^ is located near the value q — 4. The average total multiplicity is about 40. These 
values can be approximately obtained with the model parameters 

a: = 3, = 7r/12. (27) 

The only inconsistency with experiment comes from the distance between zeros of cumulant 
moments which equals 12 that is much larger than the experimentally observed ones. However, 
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it would be very naive to expect even this quite exotic quantity to be perfectly reproduced by 
such a simple phenomenological model. Further refinements of the expression (13) are possible 
to get better agreement (e.g., the nonlinear dependence of the cosine argument on r). 

Another problem appears for hadronic reactions if the same parameters are used. Particle 
production is often interpreted in terms of reggeon exchanges. In [8] it is proposed to identify 
cut disconnected vacuum-vacuum diagrams as cut reggeons. Then the average number of these 
diagrams (12) for the above values of x and <p is very large (n cut ) 14 compared to the typical 
values in widely used reggeon models, and the average number of particles in a single cut 
subgraph is quite small (n)i ~ 3. Therefore, the simple identification of cut subdiagrams with 
cut reggeons looks rather improbable. 

One of the most important features of the model is the violent behavior of the generating 
function in the complex w-plane. Its extrema are positioned at 

u k = j . / 28 
2x sin 

(k is an integer number) with maxima and minima replacing one another at each subsequent k 
and separated by the distance 

Au = — ?— r. (29) 
xsm<p 

It is especially interesting that the first extremum u\ is positioned at u± — 2 for the above 
chosen values of the parameters x and <fi. The generating function reaches extremely high value 
at this point 

G(ui) = G{2) « exp[^e 6 ] > 1. (30) 
At the next minimum w 2 — 6 the generating function is very small 

G(u 2 ) = G(6) » exp[-^e 18 ] « 1. (31) 

Maxima and minima are separated by Au = 4. 

In some way the generating function reminds the grand canonical function of the statistical 
mechanics with z = u — 1 interpreted as fugacity. It was shown [21, 22] that the singularity 
of the grand canonical function at u = 2 corresponds to the phase transition in the condensed 
matter. The obtained large value of G(2) indicates that the model probably describes the 
processes with drastically varying characteristics. The point u = 2 is the accumulating point 
of zeros of the cut generating functions 

G c (u) = £ P n u n (32) 

n=l 

for n max — > oo. These zeros are located [23, 24] near the circle \u — 1| = 1 in the complex 
w-plane. Actually, only the cut generating functions are available from experiment because of 
finite n max measurable at finite energies. It is the polynomial of u normalized at the center of 
the circle \u — 1| = 1 by the condition G c (l) = 1 with zero at u — and other n max — 1 complex 
conjugated zeros near this circle. For example, for G c {n max = 3) these two zeros are at 



jap 3 p 1 - Pi 

" 2 ' s = ~2K ± 2A < 33 » 
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with 




2P 3 



P2 



< 0. 



(34) 



At larger n max , R- e U nmax increases and tends to 2 with lmu Umax coming closer to the circle 
\u — 1\ = 1. It is remarkable that accumulation of zeros is related to the very high value of 
G(2) which reminds the singular behavior at phase transitions. 

4 Conclusions 

The simple analytical model of particle production by strong external sources is proposed. All 
the moments of multiplicity distribution and its generating function are calculable analytically. 
It rarely happens that one gets them altogether. This is the priviledge of simplest distributions 
of the probability theory but not of QCD. 

The model describes qualitatively the multiplicity distribution of produced particles fitting 
its mean value and higher order moments to experimental data and QCD findings. In particular, 
the oscillations of cumulant moments, first predicted by QCD, are reproduced. The energy 
behavior of these characteristics is related to the energy dependence of the model parameters. 
Asymptotic relations help reveal the tendencies in this dependence. 

The generating function of the distribution sharply depends on the auxiliary variable u. 
This behavior is reminiscent of the grand canonical function in statistical mechanics near the 
phase transitions. The exponential function whose exponent is an exponential function by itself 
with an oscillating factor in front of it is an unique guest in physics studies. It has no analogues 
in the well known distributions of the probability theory either. 

Even though the results rely strongly on the expression (8) obtained for 3 -theory with 
strong external sources, they can be of more general validity. The ansatz (13) combined with (8) 
reproduces very tiny features of multiplicity distributions learned from QCD and experimental 
studies. It can help get some insight for developing approximate methods of computing these 
probabilities in the field theory approach. 
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